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Abstract 

We study charge transport of quantum critical points described by conformal field theories in 2+1 
spacctimc dimensions. The transport is described by an effective field theory on an asymptotically anti- 
de Sitter spacetime, expanded to fourth order in spatial and temporal gradients. The presence of a horizon 
at non-zero temperatures implies that this theory has quasi-normal modes with complex frequencies. The 
quasi-normal modes determine the poles and zeros of the conductivity in the complex frequency plane, 
and so fully determine its behavior on the real frequency axis, at frequencies both smaller and larger 
than the absolute temperature. We describe the role of particle-vortex or S-duality on the conductivity, 
specifically how it maps poles to zeros and vice versa. These analyses motivate two sum rules obeyed by 
the quantum critical conductivity: the holographic computations are the first to satisfy both sum rules, 
while earlier Boltzmann-theory computations satisfy only one of them. Finally, we compare our results 
with the analytic structure of the 0{N) model in the large- A/" limit, and other CFTs. 
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I. INTRODUCTION 

The dynamics of quantum criticalitj^ has long been a central subject in the study of correlated 
quantum materials. Two prominent examples of recent experiments are: (z) the observation of 
criticality in the penetration depth of a high temperature superconductor at the quantum critical 
point of the onset of spin density wave order^ and [ii) the criticality of longitudinal "Higgs" 
excitations near the superfiuid-insulator transition of ultracold bosons in a two-dimensional lattice.l^ 

A complete and intuitive description of the low temperature dynamics of non-critical systems 
is usually provided by their quasiparticle excitations. The quasiparticles are long-hved excita- 
tions which describe all low-lying states, and their collective dynamics is efficiently captured by a 
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quantum Boltzmann equation (or its generalizations). The Boltzmann equation then be used to 
describe a variety of equihbrium properties, such as the electrical conductivity, thermal transport, 
and thermoelectric effects. Moreover, such a method can also address non-equilibrium dynamics, 
including the approach to thermal equilibrium of an out-of-equilibrium initial state. 

A key property of strongly-interacting quantum critical systems is the absence of well-defined 
quasiparticle excitations. The long lifetimes of quasiparticles is ultimately the justification of the 
Boltzmann equation, so a priori it appears that we cannot apply this long-established method to 
such quantum critical points. However, there is a regime where, in a sense, the breakdown of 
quasiparticle excitations is weak: this is the limit where the anomalous exponent, usually called 
T], of a particle-creation operator, 0, is small (strictly speaking, creates particles away from the 
quantum critical point). The spectral weight of the (f) Green's function is a power-law continuum, 
but in the limit — )• 0, it reduces to a quasiparticle delta function. By expanding away from the 
1] limit, one can extend to the Boltzmann method to quantum critical points, and such a 
method has been the focus of numerous studies.'^H^I 

A typical example of such Boltzmann studies is the theory of transport at the quantum critical 
point of the A^-component 0^ field theory with 0{N) symmetry in 2+1 dimensions; the N = 2 
case describes the superfluid-insulator transition of Ref. iS] Conformal symmetry emerges at the 
quantum critical point and the corresponding conformal field theory (CFT) admits a finite d.c. 
charge conductivity even in the absence of translation-symmetry breaking perturbation^ (such 
as disorder or Umklapp scattering). CFTs in 2+1 dimensions have this unique property due to 
the presence of independent but symmetry-related positive and negative charge excitations. The 
Boltzmann analysis of transport was applied in the large A^ limit and the structure of the 



frequency-dependence of the conductivity, a{u), is illustrated in Fig. 1(b) The low frequency 



behavior is as expected for weakly interacting quasiparticles: there is a Drude peak whose height 
diverges as ~ A^, and whose width vanishes as 1/A^, while preserving the total weight as A^ — t- oo. 
It is not at all clear whether such a description of the low frequency transport is appropriate for 
the A^ = 2 of experimental interest: while it is true that the anomalous exponent t] remains small 
even at A^ = 2, it is definitely not the case that the thermal excitations of the quantum critical 
point interact weakly with each other. At high frequencies, u ^ T (T is the temperature), the 
predictions of the large A^ expansion for (j{uj) seem more reliable: the result asymptotes to a non- 
zero universal constant (Too whose value can be systematically computed order-by-order in the 1/A^ 
expansion, without using the Boltzmann equation. 

In this paper, we argue for a different physical paradigm as a description of low frequency 
transport near quantum critical points, replacing the quasiparticle-based intuition of the Boltz- 
mann equation. We use the description of quantum-critical transport based on the AdS/CFT 
correspondence^ to emphasize the physical importance of "quasi-normal modes" in the charge 
response function. Formally, the quasi-normal mode frequencies are the locations of poles in the 
conductivity in the lower-half complex frequency plane i.e. the poles obtained by analytically 
continuing the retarded response function from the upper-half plane (UHP) to the second Rie- 
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FIG. 1. (a) Perspective on approaches to the charge transport properties of strongly interacting CFTs 
in 2+1 dimension. The quantum Boltzmann approach apphes to the expansion of the 0{N) model: 
its starting point assumes the existence of weakly interacting quasiparticles, whose collisions control the 
transport properties. In the present paper we start from the "nearly perfect" quantum liquid obtained 
in the A'c — s- oo limit of a SU(A''c) super Yang-Mills theory, which has no quasiparticle description. 
Holographic methods then allow expansion away from this liquid (A is the 't Hooft coupling of the gauge 
theory), (b) Structure of the charge conductivity in the quantum Boltzmann approach. The dashed line 
is the N = oo result: it has a delta function at zero frequency, and a gap below a threshold frequency. 
The full line shows the changes from corrections, (c) Structure of the charge conductivity in the 
holographic approach. The Nc = oo result is the dashed line, and this is frequency independent. The full 
line is the conductivity obtained by including four-derivative terms in the effective holographic theory for 
7 > 0. 



mann sheet in the lower-half plane (LHP). By considering a particle- vortex dual (or "S-dual") 
theory whose conductivity is the inverse of the conductivity of the direct theory, we also associate 
quasi-normal modes with the poles of the dual theory, which are the zeros of the direct theory. 
Both the pole and zero quasi-normal modes are directly accessible in AdS / CFT methods,'^^^ and 
are related to the normal modes of excitations in the holographic space: the normal modes have 
complex frequencies because of the presence of the "leaky" horizon of a black brane; see Fig. |2} 
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FIG. 2. AdS spacetime with a planar black brane. The current ( J^) correlators of the CFT are related 
to those of the U(l) gauge field (A^) in the AdS (bulk) spacetime. The temperature of the horizon of 
the black brane is equal to the temperature of the CFT. The horizon acts as a "leaky" boundary to the 
bulk normal modes, which consequently become quasi-normal modes with complex frequencies. These 
quasi- normal modes specify the finite temperature dynamic properties of the CFT. 

We will show that a knowledge of these modes allows a complete reconstruction of the frequency 
dependence of the conductivity, (t{u), extending from the hydrodynamic regime with u <^ T, to 
the quantum critical regime with u ^ T. Moreover, these quasi-normal mode frequencies are 
also expected to characterize other dynamic properties of the quantum critical system: the recent 
work of Bhaseen et alJ^ showed that the important qualitative features of the approach to thermal 
equilibrium from an out-of-equilibrium thermal state could be well understood by a knowledge of 
the structure of the quasi-normal mode frequencies. 

Apart from the quasi-normal modes, the long time dynamics also exhibits the well-knowrP^ 
classical hydrodynamic feature of 'long time tails' (LTT). The LTT follow from the principles 
of classical hydrodynamics: arbitrary long wavelength hydrodynamic fluctuations lead to the al- 
gebraic temporal decay of conserved currents. The LTT depend only upon various transport 
coefficients, thermodynamic parameters, and a high frequency cutoff above which hydrodynamics 
does not apply. In the quantum-critical systems of interest here, this high frequency cutoff is 
provided by the quasi-normal modes. Thus the LTT describe the dynamics for frequencies a; -C T, 
while the quasi-normal modes appear at w ~ T and higher. We emphasize that the value of the 
d.c. conductivity, a{uj/T = 0), is determined by the full CFT. The non-analytic small frequency 
dependence associated with the LTT can be obtained from the effective classical hydrodynamic 
description which takes the transport coefficients of the CFT treatment as an input. The focus of 
the present paper will be on the quasi-normal modes, and we will not have any new results on the 
LTT; the description of the LTT by holographic methods requires loop corrections to the gravity 
theory,!^ which we will not consider here. 
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From our quasi-normal mode perspective, we will find two exact sum rules that are obeyed by 
the universal quantum critical conductivity, cr{uj), of all CFTs in 2+1 dimensions with a conserved 
U(l) charge. These are 

/■oo 

du[^a{u) - a^] = 0, (1) 
1 1 







oo 



duj 



3f? 



0. (2) 



Here Uoo is the limiting value of the conductivity for a; ^ T (in applications to the lattice models 
to condensed matter physics, we assume that u always remains smaller than ultraviolet energy 
scales set by the lattice). The first of these sum rules was noted in Ref. |22l From the point of 
view of the boundary CFT, Eq. Q is quite natural in a Boltzmann approach: it is similar to 
the standard /-sum rule, which we extend to CFTs in Appendix [A} There we connect it to an 
equal-time current correlator, which we argue does not depend on IR perturbations such as the 
temperature or chemical potential. The second sum rule follows from the existence of a S-dual 
(or "particle- vortex" dual) theory^^ 18 1 231 B5] -^j^ose conductivity is the inverse of the conductivity of 
the direct theory. Although it can be justified using the direct sum rule, Eq. ([T]), applied to the S- 
dual CFT, whose holographic description in general differs from the original theory, we emphasize 
that it imposes a further constraint on the original conductivity. To our knowledge, the second 
sum rule has not been discussed previously. All our holographic results here satisfy these two 
sum rules. We show in Appendix [B| that the N = 00 result of the 0{N) model in Ref. |3| obeys 
the sum rule in Eq. ([T|, a feature that was not noticed previously. However such quasiparticle- 
Boltzmann computations do not obey the sum rule in Eq. ([2]). The holographic computations of 
the conductivity are the first results which obey not only the sum rule in Eq. ([T|, but also the 
dual sum rule in Eq. ([2]). 

In principle, the quasi-normal mode frequencies can also be determined by the traditional meth- 
ods of condensed matter physics. However, they are difficult to access by perturbative methods, or 
by numerical methods such as dynamical mean-field theory.'^SI One quasi-normal mode is, however, 
very familiar: the Drude peak of quasiparticle Boltzmann transport, appearing from the behavior 
a{uj) ~ cro/(l ^ ioor), corresponds to a quasi-normal mode at u = —i/r. In a strongly-interacting 
quantum critical system, we can expect from the arguments of Ref. |4] that this peak would trans- 
late to a quasi-normal mode at a; ~ —iT. As we will see in detail below, this single Drude-like 
quasi-normal mode does not, by itself, provide a satisfactory description of transport, and we 
need to understand the structure of the complete spectrum of quasi-normal modes. And the most 
convenient method for determining this complete spectrum is the AdS/CFT correspondence. 



As we indicate schematically in Fig. 1(a), the AdS/CFT description becomes exact for certain 
supersymmetric gauge theories in the limit of a large number of colors, Nc, in the gauge group.'^^'^ 
This theory has no quasiparticles, and in the strict iV^ = 00 limit the conductivity is frequency 



independent even at T > 0, as indicated in Fig. 1(c) Our quasi-normal mode theory expands 
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away from this frequency-independent limit, in contrast to the free particle limit of the Boltzmann 
theory (in the latter limit, the Drude contribution becomes cr{u) ~ T6{uj)). We describe the basic 
features of o"(c<j) obtained in this manner in the following subsection. Because strong interactions 
are crucial to the structure of <j{ou) at all stages, and there is no assumption about the existence of 
quasiparticles, we expect our results to be general description of a wide class of strongly interacting 
quantum critical points. 



A. Generic features of the finite-T conductivity of a CFT 

The frequency dependent conductivity of a CFT in 2+1 dimensions at finite temperature will 
naturally be a function of the ratio of the frequency to the temperature, u/T, which we will denote 
as w, with a factor of An convenient in the holographic discussion, 

^ = 4^- (3) 

In general, we do not expect the conductivity of a generic CFT to be a meromorphic function 
of the complex frequency w, i.e. analytic except possibly at a discrete set of points where it 
has finite-order poles, all in the LHP. (The latter condition follows from the causal nature of the 
retarded current-current correlation function.) The absence of meromorphicity for the conductivity 
of an interacting CFT, or the presence of branch cuts, can be attributed to the LTT.™Ql the 
present paper we will not discuss LTT and focus on the meromorphic structure of the conductivity. 
On the one hand, such a description should be valid for CFTs that have a holographic classical 
gravity description.!^ For example, there is strong evidence that certain super Yang-Mills large- A*";, 
gauge theories are holographically dual to classical (super)gravity and do not have LTT, which 
are suppressed by l/A^^ compared to the leading meromorphic dependence On the other hand, 
we believe that understanding the meromorphic structure is a first step to understanding the full 
analytic structure of generic CFTs, and do not expect branch cuts from the LTT to significantly 
modify the poles and zeros of the quasi-normal modes at frequencies of order T or larger. 

The meromorphic condition is tantamount to assuming that in response to a small perturbation, 
the system will relax exponentially fast to equilibrium at finite temperature. In addition to LTT, 
we expect deviations from such behavior to occur at a thermal phase transition for instance, where 
power law relaxation will occur. In that case a is not expected to be meromorphic and branch 
cuts can appear. Another exception is free CFTs, such as the 0{N) model in the limit where 
N ^ oo, where we find poles and zeros directly on the real frequency axis, as well as branch cuts, 
as shown in Fig. 3(c)[ We restrict ourselves to the finite temperature regime of an interacting 



conformal quantum critical point with a classical gravity description and do not foresee deviations 
from meromorphicity 
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Moreover, we expect the universal conductivity to go to a constant as w ^ cxoPEH 



a{w — 7- oo) = (Too < oo , w E M. . (4) 

Such a well-defined limit will generally not exist as one approaches complex infinity along certain 
directions in the LHP. This is tied to the fact that a will not necessarily satisfy the stronger 
condition of being additionally meromorphic at infinity. In other words, s{z) := a{l/z) is not 
necessarily meromorphic in the vicinity of the origin, z = 0. If it were, cr{w) would be a rational 
function, the ratio of two finite-order polynomials, and would have a finite number of poles (and 
zeros). In our analysis, we shall encounter a class of CFTs whose conductivity has an infinite set 
of simple poles, and is thus not meromorphic on the Riemann sphere CU {oo}. A familiar example 
of such a function is the Bose-Einstein distribution, nsiw) = l/{e^ — 1), which is meromorphic, 
but not at infinity because it has a countably infinite set of poles on the imaginary axis. In fact 
nsil/z) has an essential singularity at 2; = 0. 

A further generic property that a satisfies in time-reversal invariant systems is refiection sym- 
metry about the imaginary-frequency axis: a{—w*) = (7{w)*, which reduces to evenness or oddness 
for the real and imaginary parts of the conductivity at real frequencies, respectively. In particular, 
this means that all the poles and zeros of a either come in pairs or else lie on the imaginary axis. 
Following this discussion, we can express the conductivity as 

^ ' n poles U,{^ - UJ^ - ^m){w + vr;;) ' ^> 

where ^ denotes zeros and vr poles; {C/*, vr^} and {Cm, T^m} he on and off the imaginary axis, respec- 
tively. In this sense the poles and zeros contain the essential data of the conductivity. Actually, 
since a{w — oo)/croo = 1 on the real axis, which also holds for all directions in the UHP, they 
entirely determine o joao- In the current holographic analysis, all the poles and zeros are simple, 
excluding double and higher order poles. We suspect this is a general feature of correlated CFTs. If 
one is interested in the behavior on the real frequency axis only, the expression for the conductivity 
arising from the AdS/CFT correspondence can be truncated to a finite number of poles and zeros: 
we will show in Section |TTE that this leads to reasonable approximations to the conductivity on the 



real frequency axis. Such a truncated form can be compared with experimentally or numerically 
measured conductivities for systems described by a conformal quantum critical point. 

As we will show in this paper, the holographic methods allow easy determination of the poles 
in the conductivity, which are identified as the frequencies of the quasi-normal modes of the 
theory on AdS4 in the presence of a horizon at a temperature T. Moreover, the zeros in the 
conductivity emerge as the frequencies of the quasi-normal modes of a S'-dual (or "particle- vortex" 
dual) theory.'^S'ISESCSl '\Ye summarize our holographic results for a particular parameter value in 
Fig. |3| along with the corresponding results for the 0{N) model at = 00. The 0{N) model has 
a pole at w = 0, corresponding to the absence of collisions in this model at A^ = 00. This turns 
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FIG. 3. (a) Poles (crosses) and zeros (circles) of the holographic conductivity at 7 = 1/12. (b) Real and 
imaginary parts of the holographic conductivity on the real frequency axis, (c) Poles and zeros of the 
0{N) model at = 00; the zeros coincide with branch points, and the associated branch cuts have been 
chosen suggestively, indicating that the branch cuts transform into lines of poles and zeros after collisions 
have been included, (d) Conductivity of the 0{N) model at = 00; note the delta function in the real 
part at w = 0, and the co- incident zero in both the real and imaginary parts at a; = 2 A. In these figures 
A/T = 21n((-v/5 + l)/2), and the 0{N) computation is reviewed in Appendix |b[ 



into a Dru de-like pole on the imaginary axis, closest to the real axis in the holographic result. We 
show in Appendix |B] that the 0{N) model also has a pair of zeros on the real axis, and this is 
seen to correspond to zeros just below the real axis in the holographic result. Finally, the 0{N) 
model has a pair of branch points on the real axis; the location of the branch cuts emerging from 
these branch points depends on the path of analytic continuation from the upper half plane. We 



have chosen these branch cuts in a suggestive manner in Fig. 3(c), so that they correspond to the 
lines of poles and zeros in the lower-half plane of the holographic result. So we see a natural and 
satisfactory evolution from the analytic structure of the collisionless quasiparticles of the 0{N) 
model, to the quasi-normal modes of the strongly interacting holographic model. 

The outline of our paper is as follows. The holographic theory on AdS4 will be presented 
in Section [TT| We will use the effective field theory for charge transport introduced in Ref. EH 
expanded to include terms with up to 4 spacetime derivatives. The quasi-normal modes will be 
computed using methods in the literature J ^^ * ^^ * ^^ * ^ Section III will turn to the traditional quantum 
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Boltzmann methods where new results regarding the analytic structure are given; in particular, 
we find that the low frequency Boltzmann conductivity can be accurately represented by a single 
Drude pole. 



II. HOLOGRAPHIC ANALYSIS 

The AdS/CFT holographic correspondence we use arose from the study of non-abelian su- 
persymmetric gauge theories in the limit of a large number of colors, for example with gauge 
group SU(Ai'c), Nc — oo. By taking an appropriate limit for the gauge coupling, such theories 
are strongly interacting yet they can be described by weakly coupled gravity in an Anti-de-Sitter 
(AdS) spacetime with one extended additional spatial dimension, and six or seven compactified 
ones. The fixed-point CFT describing the strongly correlated gauge theory can be seen as existing 
on the boundary of AdS. Different correlation functions on the boundary quantum CFT, such as 
the charge- current ones of interest to this work, can be computed by using the bulk (semi-) classical 
gravitational theory. For instance, the current operator corresponding to a global U(l) charge in 
the CFT can be identified with a U(l) gauge field in the higher dimensional gravitational bulk 
(Fig. |2|. We refer the reader to a number of review^* ^^*^^ ! with condensed matter applications in 
mind and proceed to the holographic description of transport in 2+1 dimensional CFTs. 

These CFTs are effectively described by a gravitational bulk theory in 3+1 dimensions. In the 
case of the supersymmetric AB JM modeP^ in a certain limit with an infinite number of colors, the 
holographic dual is simply Einstein's general relativity in the presence of a negative cosmological 
constant resulting in an AdS4 spacetime. Charge-transport correlations functions in the CFT can 
be obtained from those a U(l) probe gauge field with Maxwellian action in the AdS background. 
It was showEp' that the conductivity of the large- A^c ABJM model is frequency independent due 
to an emergent S-duality. Ref. EH discovered that deviations from self-duality are obtained by 
considering 4-derivative corrections to the Einstein-Maxwell theory, which can potentially arise at 
order 1/A in the inverse 't Hooft coupling. The effective action for the bulk gravitational theory 
discussed in Ref. EU reads 



5'buik = / d'^Xy/^ 



(6) 



where g is the determinant of the metric Qab with Ricci scalar R; F""^ is the field strength tensor 
of the probe U(l) gauge field Aa holographically dual to the current operator of a global charge 
of the CFT. (We use Roman indices for the 3+1 spacetime, and Greek ones for the boundary 
2+1 spacetime.) Such an action was also considered in Ref. ESI The 4-derivative contribution 
to charge-transport can be encoded in the last term, proportional to 7. Cabcd is the (conformal) 
Weyl curvature tensor; it is the traceless part of the full Riemann curvature tensor, Rabcd'- Cabcd = 
Rabcd-{ga[cRd]b- gb[cRd]a) + lRga[cg d]b- We observe that the 7-term directly couples the probe U(l) 
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gauge field to the metric. L is the radius of curvature of the AdS4 space while the gravitational 
constant is related to the coefficient of the two-point correlator of the stress-energy tensor T^y 
of the boundary CFT (for a review, see Ref. I39ll . an analog of the central charge of CFTs in 
1-l-lD. The gauge coupling constant gl = l/o"oo dictates the infinite-w conductivity, which we 
shall set to 1 throughout, effectively dealing with ajooo- The crucial coupling in this theory is 
the dimensionless parameter controlling the four-derivative term, 7; it determines the structure 
of a three-point correlator between the stress-energy tensor and the conserved current. Stability 
constraints in the theory imp Ij^^ that I7I < 1/12, and we explore the full range of allowed 7 values 
here. Positive values of 7 yield a low-frequency peak in the conductivity as shown in Fig. |l(c) or 



Fig. 7(a) , while negative values of 7 give rise to a low-frequency dip illustrated in Fig. |7(b) 



as 



may be expected from a theory of weakly interacting vortices. Explicit computations of 7 directly 
from the CFT yield valued^ in line with these expectations. 



In the spirit of the effective field theory approach of Ref. EH we should also consider adding 
other terms to Eq. ([6]) involving fields other than Fah and the metric tensor .SSI The most important 
of these are possible "mass" terms which tune the CFT away from the critical point at T = 0. 
Such terms are not present in the CFT at T = 0, but their values at non-zero T are precisely such 
that the expectation value of the mass operator does not change: e.g. in the quantum critical 
0{N) model of Appendix |a| (0^) is T- independent.'^ The mass terms can be included in the 
holographic theory by allowing for a scalar dilaton field, $, and this can modify charge transport 
via a term ~ (^F^yF^^ . In the holographic theory, in the absence of external sources, such a dilaton 
does not acquire an expectation value at T > when it is not present at T = 0. And external 
sources coupling to the gauge field only modify $ at quadratic order, and so $ can be neglected 
in the tree-level linear response. Thus even after allowing for additional fields, 7 remains the only 
important coupling determining the structure of the charge transport at non-zero temperatures. 



In the absence of the gauge field, which is here only a probe field used to calculate the linear 
response, the metric that solves the equation of motion associated with Sbuik is: 

ds' = T2{-mde + dx- + dy') + ^y (7) 

where /(r) = 1 — rg/r^, and r is the coordinate associated with the extra dimension. The CFT 
exists on the boundary of AdS, r — 00, on the Minkowski spacetime parameterized by {t,x,y). 
We emphasize here that the holographic theory is naturally written in real time allowing direct 
extraction of the retarded current-current correlation function characterizing the conductivity. 
Eq. ([7]) corresponds to a 3+lD spacetime with a planar black hole (BH) whose event horizon 
is located at r = tq, and that asymptotically tends to AdS4 as r — )■ 00. We thus refer to it 
as Schwarzchild-AdS, or S-AdS. The position of the event horizon is directly proportional to the 
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temperature of the boundary CFT, 



As T — )• 0, the black hole disappears and we are left with a pure AdS spacetime, which is holo- 
graphically dual to the vacuum of the CFT. The statement that the thermal states of the CFT 
can be accessed by considering a BH in AdS can be heuristically understood from the fact that 
the BH will Hawking radiate energy that will propagate to the boundary and heat it up. 

It will be more convenient to use the dimensionless coordinate u = ro/r, such that Eq. ([T]) 
becomes 

ds' = {-f{u)dt' + dx' + dy') + , f{u) = l-u\ (9) 

The boundary, r = cxd, is now at u = 0, while the BH horizon is at u = 1. 

The equation of motion (EoM) for the probe gauge field is the modified Maxwell equation 

Va(F'^''-47L2C"^'"='^F,,) = 0, (10) 

where Va denotes a covariant derivative with respect to the background metric, Qab- As we are 
interested in the current correlator in frequency-momentum space, we Fourier transform the gauge 
field: 

/d-^k 
j^,e-^-'^^'^-^Aaico,K,ky,u), (11) 

where the coordinate u was left un-transformed since there is no translational invariance in that 
direction. We shall actually solve for the full u-dependence of Aa. We work in the radial gauge 
Au = 0. Without loss of generality, we also set the spatial momentum to be along the x-direction, 
{kx,ky) = {k,0). In the limit where /c — )■ 0, appropriate to a uniform "electric" field coupling to 
the global charge, the equation of motion for the transverse component. Ay, reads 

h' 9w'^ 

K + J^K + jr^y = ^^ (12) 

where we have defined the dimensionless frequency w in Eq. (|3|), and primes denote derivatives 
with respect to u. The function h{u) is simply fg, where g = 1 + 4'~fu^ takes the same form as 
f = 1 — u^. As g{u) fully encodes the 7-dependence, we wish to make its role more transparent 
by rewriting the above equation: 



Ay+(^ + '-)A'y+'-^Ay = 0. (13) 
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The term g' /g = 127u^/(l + 4'~fu^) is seen to be proportional to 7, and as such, goes to zero as 
u — > consistent with the fact that the Weyl tensor vanishes in the pure AdS spacetime, which is 
said to be conformally-flat. 

The AdS/CFT correspondence provides an expression for the conductivity of the CFT in terms 
of the transverse gauge field auto-correlator evaluated at the boundary, u = 0, 



yy 



UJ 



(14) 



n=0 



where cr(a;) is the complex valued conductivity, and Qyy{u,u) is the retarded Ay auto-correlation 
function. More specifically one 

gets™ 



i d^Ay 



3w A,, 



(15) 



u=0 



where Ay solves the equation of motion Eq. (13) with suitable boundary conditions, as discussed 



below. The above equation, central to our analysis, has the following heuristic explanation: Ay{0) 
acts as a source for the current, while duAy{0) is the corresponding response. We will see in 



Section [II C| that the quasi-normal modes, i.e. the poles of conductivity in the LHP, correspond 
to driving frequencies at which a "response" exists in the limit of vanishing source strength. 



A. Direct solution of conductivity 



The real part of the conductivity on the real frequency axis (retarded correlator) was numerically 
obtained in Ref. |211 We extend their analysis from real to complex frequencies, w E C The 
boundary conditions necessary to solve Eq. (13) are imposed at the BH event horizorP^^at u = 1. 



To obtain them we examine the EoM near the horizon, which admits the following two solutions: 
These correspond to outgoing and ingoing waves from the point of view of 



Ay ~ 



1 



u 



the BH, respectively. The retarded correlator is obtained by choosing the ingoing condition. To 



implement this in the numerical solution, we factor out the singular behavior: Ay 



-u 



'F{u) 



where F{u) is the sought-after function; it is regular at the horizon. From Eq. (15), we see that 



we are free to fix one of the two boundary conditions, either for Ay{l) or ^4^(1), to an arbitrary 
finite constant without altering the conductivity. We impose Ay{l) = F{1) = 1. The appropriate 
boundary condition for F' can be obtained by examining the differential equation near m = 1 as is 
discussed in Ref. EU and in Appendix [C| 

All the poles of the conductivity are in the LHP, as it is obtained from the retarded current- 
current correlation function. The numerical result is shown in Fig. 4(a) and Fig. 4(c)| for the 
two values of 7 saturating the stability bound, 7 = ±1/12, respectively. Fig. 4(a) shows the 
conductivity for 7 = 1/12, which corresponds to particle-like transport with a Drude peak at small 
real frequencies as can be seen on the real w-axis, or more clearly in Fig. 3(b) or Fig. 7(a) Such 
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(a)3?{CT(w;7 = 1/12)} 



(b)3?{o-(w;7 = 1/12)} 




(c)3fi{cr(w;7 = -1/12)} {d)n{a{w;j = -1/12)} 



FIG. 4. Conductivity a and its S-dual a = 1/a in the LHP, w" = < 0, for I7I = 1/12. The zeros 
of a{w,'y) are the poles of d-{w;^). We further note the quahtative correspondence between the poles of 
a^Wjj) and the zeros of a{w; —7). 



low-frequency behavior is dictated by a Drude pole, located closest to the origin. The numerical 
solution also shows the presence of satellite poles, the two dominant ones being shown. These 
are symmetrically distributed about the '^w axis as required by time-reversal, and are essential 
to capture the behavior of cr beyond the small frequency limit. In contrast, the conductivity at 
7 = —1/12 in Fig. 4(c) shows a minimum at w = on the real axis, see also Fig. 7(b) for a plot 
restricted to real frequencies. The corresponding pole structure shows no poles on the imaginary 
axis, in particular no Drude pole. The conductivity at 7 = —1/12 is said to be vortex-like because 
it can be put in correspondence with the conductivity of the CFT S-dual to the one with 7 = 1/12, 
as we now explain. 



B. S-duality and conductivity zeros 



Great insight into the behavior of the conductivity can be gained by means of S-duality, a 
generalization of the familiar particle- vortex duality of the 0(2) model. S-duality on the boundary 
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CFT is mirrored by electric-magnetic (EM) duality for the bulk U(l) gauge field, which we now 
briefly review. Given the abelian gauge theory for the U(l) bulk field Aa, we can always perform 
a change of functional variables in the partition function to a new gauge field Aa by adding the 
following term to S'buik) Eq- (js]): 

S' = J (fx^g ^e'^'^'^Aad^F,, , (16) 

with the corresponding functional integral for Aa- Performing the integral over Aa would simply 
enforce the Bianchi identity, e°'^^'^di,Fcd = 0, implying Fab = daA^ — di,Aa, where Sabcd is the fully- 
antisymmetric tensor in 3-l-lD with Stxyu = \f-~9- 

If instead one integrates out Aa first, a new action in terms of Aa results: 

^buik = - / d'x^g -^FabX^'^'^F,, , (17) 



89i 

where we have defined the field strength of the dual gauge field. Fab = daAb — dbAa, and dual 
coupling g4 = l/g^. An exactly analogous action holds for Aa without the hats. The rank-4 
tensors X, X are shorthands to simplify the actions: 

Xab'^ = Iab"'-8lL^Cab''\ (18) 

Xab'^ = \eab'^{X-%f^egh'\ (19) 

with the rank-4 tensor lab'^'^ = Sa'^Sb'^ — Sa'^Sb'^, the identity on the space of two-forms, for e.g. 
^ab = ^lab'^'^'Fcd- The iuvcrsc tensor of X is then defined via ^{X'^^)ab'^'^Xc/-l' = lab'^^ ■ In terms of 
the X-tensors, the EoM for Aa and Aa simply read: 

Vfe(X'^^^'^F,,) = 0, (20) 
Vfe(X'^'^'^F,,) = 0. (21) 

It can be shown^^ that for small 7, the dual X-tensor has the following Taylor expansion: 

Xab"" = lab"" + 87^2^^'='^ + 0(7') , (22) 

= Xab''" |,^_, + 0(7') • (23) 

We thus see that if 7 = 0, X = X and the actions, and associated EoM, for A and A have the 
same form. In that case, the two theories are related by an exchange between electric and magnetic 
fields: the standard EM (hodge) self-duality of electromagnetism. In contrast, in the presence of 
the 4-derivative term parameterized by 7, the EM self-duality is lost. However, at small 7 the EM 
duality is particularly simple and will serve as a guide for any finite 7: the holographic theory for 
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7 maps to the one for —7, neglecting 0(7^) contributions. 

Let us now examine the impact of this bulk EM duality, A A, on the boundary CFT. The 
holographic correspondence relates the bulk gauge field A to the current of a global U(l) charge 
of the CFT, J. In the same way, the dual gauge field A will couple to the current J of the S-dual 
CFT, which generically differs from the original CFT. Just as the conductivity of the original CFT, 
a, is related to the J auto-correlator, the conductivity of the S-dual CFT, a, will be obtained from 
the J auto-correlator. The conductivities of the S-dual CFT pair are in fact the inverse of each 
other: 

a{w;^) = -^ (24) 
cr(w;7) 

where we emphasize that this relation holds for the complex conductivities, a = 3fJcr -|- iQa. We 
present the short proof here using results of Ref. IM1 (We note that such a result was derived 
for a specific class of CFTs in Ref. [TSi) We begin with the general form of the retarded current- 
current correlation function: Q^i^{u, q) = \/ q\q^{P'^^K'^ ibj , q) + Pl^^K^ {uj , q)), with the orthogonal 
transverse and longitudinal projectors P^'^: = = P^ = 0, P^ = 6ij — qiqj/q'^, and by 
orthogonality: = [r^^i, — q^qv / {qxq^)] — The Minkowski metric was introduced, rj^y = 

diag(— 1, 1, 1), such that q\q'^ = rfwiq^q^' = —ij^ -|- g^. Of interest to us is the holographic relation 
between the transverse correlator giving the conductivity and the bulk gauge field correlator, Qy,^: 

^q^ -UJ^K^{uj,q) = Gyyiu, q)\^^Q = uja{uj, q) /i , (25) 

where (j{u, q) is the frequency and momentum dependent conductivity. The same expression (with 
hats) holds in the S-dual theory. Using the action of EM duahty on the bulk, Ref. [2l] showed the 
relation: 

K^{u,q)k'^{uj,q) = l, (26) 

that relates the transverse current-current correlator of the original CFT to the longitudinal one 
of the dual CFT. When combined with the fact that in the limit of vanishing spatial momentum, 
g — )■ 0, rotational invariance enforces K'^{u, q) = K^{uj, q), which is also naturally true with hats, 
we obtain 



By virtue of Eq. (25) and its dual version, this concludes the proof of Eq. (24). 



The poles of the dual conductivity, a = 1/a, then must correspond to the zeros of the conduc- 
tivity, a, and vice versa. As a consequence, we see that S-duality interchanges the locations of the 
conductivity zeros and poles. This is is consistent with the direct solution shown in Fig. |4} Take for 
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example the theory at 7 = 1/12, Fig. |4(a) it will have a Drude pole on the imaginary axis, which 



gives rise to a Drude peak at small frequencies. Under S-duality this pole becomes a Drude zero 
of a, Fig. |4(b) , and the conductivity of the new theory will have a minimum at small frequencies. 



As we saw above, changing the sign of 7 corresponds to an approximate S-duality valid for 
I7I ^ 1. More generally, in terms of the "pole/zero-topology" or ordering, both operations are 
equivalent. Indeed, if we consider the pole/zero structure of the positive frequency branch of the 
conductivity > (which is sufficient by time-reversal) and order the poles and zeros according 
to their norm, we get the following two equivalence classes: 

pole — zero — pole — zero — . . . — t- particle -like (7 > for e.g.) , (28) 
zero — pole — zero — pole — . . . — )■ vortex-like (7 < for e.g.) , (29) 

where the first label (in bold) designates the Drude pole or zero. Both S-duality and 7 — )■ —7 
interchange these two analytic structures. This underlies the qualitative correspondence between 



the pole structure of a{w;'y) and that of a{w; —7); for example, compare Fig. 4(a) and Fig. 4(d) 



or Fig. 4(c) and Fig. 4(b) The correspondence quantitatively improves in the limit of small 7. 
Explicitly, 

a(w;7)^— ^ V' (30) 

a{w; -7) 

holds because performing a — )■ l/o" together with 7 — )■ —7 is approximately tantamount to two 
S-duality transformations and is equivalent to the identity, modulo C(7^) terms. 

Finally, we mention that for a given 7 it is not possible to find a 7' such that a{w; 7) = 
(t(w;7'). In other words, the dual of the boundary CFT with parameter 7 cannot correspond to 
the original CFT with a different parameter 7'. This can be seen as follows. We first require that 
the relation hold true at zero frequency: o"(0; 7) = cr(0; 7'), which implies 1/(1 + 47) = 1 -|- 47' or 
7' = (jTij^ — l)/4, where we have used cr(0;7) = 1 + 47 (see Refs. IMl and 1381) . Although for this 
value of 7', a{w; 7) and a{w; 7') agree for both w,l/w = 0, v/e have numerically verified that they 
always disagree at intermediate frequencies, the disagreement decreasing as 7 — )■ 0, in which limit 
7' ^ —7. The absence of a 7' satisfying a{w; 7) = cr(w; 7') is in accordance with the fact that 
holgraphic action of the S-dual CFT contains terms beyond C abcdF""^ F*^"^ . The latter is only the 
first term in the Taylor expansion in 7. 

We now turn to a better method of determining the poles and zeros, as the direct solution of 



Eq. (12) can only reliably capture the poles nearest to the origin. The main problem with the 



direct solution of the differential equation for Ay, Eq. (12), is that the Fourier modes Ay{u]w) at 



the UV boundary, m = 0, generically grow exponentially as the imaginary part of the frequency 
becomes more and more negative making the numerical results unstable. Although an exception 
occurs at the poles, where Ay{u = 0; cUpoie) vanishes (see below), it is hard to untangle the true 
analytical structure from the numerical noise, hence the need for a more sophisticated approach. 
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C. Quasi-normal modes and poles 



We present an alternative and more powerful method of capturing the poles by considering the 
so-called quasi-normal modes (QNMs) of the gauge field in the curved S-AdS4 spacetime. These 



modes are eigenfunctions of the EoM for Ay, Eq. (12): 



^'n + jA. + ^ A = , (31) 

where An is a QNM with frequency W7„. The QNM have the special property that they vanish at 



the boundary: An — )■ as m — )■ 0. From the expression for the conductivity, Eq. (15), we can see 
that this will lead to Wn being a singular point of the conductivity: 



d A 

a{wn) 



A 



^"•^"(°)^±oo. (32) 



where duAn{0) is generically finite. The name quasi-normal instead of normal is used because 
the eigenfunctions An diverge approaching the BH horizon, u = 1. This follows from the above- 
mentioned asymptotic form near the horizon. An ~ (1 ~ m)"™" = (1 — implying a 
divergence for frequencies in the LHP. As predicted by the AdS/CFT correspondence and verified 
by our numerical analysis, shown in Fig. |5j the QNMs indeed agree with the poles of the conduc- 
tivity shown in Fig. |4] and more precisely in Fig. 11 The QNMs are found by using a Frobenius 



expansion 

M 



A, = M/(M)-™^a„(M-Mr, (33) 



m=0 



where we have factored out the behavior near the event horizon, f{u) ™ ~ (1 — -u) ™, and near 
the boundary, u. We have chosen to Taylor expand around 'U = l/2;M + lis the number of 



terms in the truncated series. Substituting Eq. (33) in Eq. (12) yields a matrix equation for the 
coefficients, am'. 



M 



J2 Bimam = , (34) 



m=0 



where the l.h.s. is the coefficient of {u — u)\ < / < M. Note that Bim = Bimiw) and am = a-miw) 
both depend on the frequency, and although not explicitly shown, on 7 as well. For fixed 7, 
this homogeneous system of linear equations has a solution at a set of frequencies {wn} at which 
det B{wn) = 0. Or equivalently, when the smallest-normed eigenvalue of B, Amin, vanishes, which 
we find more convenient to implement numerically. Plots of l/|Amin| (multiplied by an exponential 
function to improve the visibility) as a function of w are given in Fig. Islfor I7I = 1/12. The QNMs 
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w' 





w' 



(a)a(7 = 1/12) 



(b)a(7=l/12) 




(c)a(7 = -1/12) 



(d)a(7 = -l/12) 



FIG. 5. Quasi-normal modes (bright spots) of the transverse gauge mode for 7 = |1/12| in the complex 
frequency plane, w = w'+iw". The QNMs correspond to the poles of the conductivity (a &: c). EM duality 
yields the QNMs of the dual gauge mode, and these correspond to the poles of the dual conductivity, 
a{w) = l/a{'w), i.e. the zeros of (7{w), see panels b & d. 



are the bright spots. In obtaining the QNMs of the dual conductivity, a = 1/a, we have used the 
EoM for the dual gauge field i, Eq. (|2T|: 

A:- 



f 9 



A' 



0. 



(35) 



It differs from the one for Ay, Eq. (13), by the negative sign.^^' Note that this shows that 7 — t- —7 
does not exactly correspond to S-duality, because the former would give —g' /{I — A'-^v?) 7^ —g'/g, 
where g = 1 + A'^u^ . 

Whereas the direct solution only gives reliable answers up to '^w ~ — 1, the QNM approach has 
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(c)7 ~ 0+ 

FIG. 6. Illustration of the motion of the poles and zeros as 7 goes to zero, in 3 steps: 7 = 10~^ — )• 10~^ — )> 
10~^ — )• O"*". In each panel the motion is from bold to thin as 7 decreases; with crosses representing poles 
while circles, zeros, a) Blue and thick markers are for 7 = 10~^, while the red an thin ones for 7 = 10~^. 
b) Red and thick markers are for 7 = 10~^, while the green an thin ones for 7 = 10~^. c) "Zipped" 
pole-zero structure for 7 ~ O"*", where only poles and zeros far from the origin will lie off the imaginary 
axis. 



a wider range of applicability and is numerically more stable giving us more insight into the analytic 
structure. We have performed a WKB analysis in Appendix |D] to determine the asymptotic QNMs 
for \w\ ^ 1. We next examine the transition that occurs when going from positive to negative 
values of 7. 



D. Pole motion and S-duality 

The motion of the poles and zeros as 7 changes sign is illustrated in Fig. [6] for 7 > 0. For 7 < 0, 
one simply interchanges the zeros and poles, i.e. the crosses and circles. The pole/zero motion can 
be loosely compared with a "zipper mechanism". The arrows in Fig. |6] show the non-trivial motion 
of a pair of poles or zeros as they become "zipped" to the imaginary axis. (A caveat regarding the 
arrows: by time-reversal symmetry, w — )■ — w*, so we cannot say which pole goes to which once 
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they become pinned to the imaginary axis. The arrows are just a guide.) For sufficiently small 7, 
each point on the imaginary axis located at w^^ = —in/2, where n is a positive integer, will have 
a pole and zero arbitrarily close to it. When 7 = 0, they will "annihilate" as it should because the 
complex conductivity for 7 = has no poles or zeros as it takes the constant self-dual value for all 
complex frequencies. It should be noted that since w = u/AirT, the annihilation frequencies are 

uj^^P = -i27mT , n = 1, 2, 3, . . . , (36) 

i.e. the bosonic Matsubara frequencies in the LHP. Although this results seems natural, we do not 
have a clear explanation for it and leave the question for future investigation. Finally, from the 
direct numerical solution of the EoM, we have looked at the residue of the pole near w = —i/2 
(closest to the origin), and have found that it decreases linearly with 7, consistent with the 7 = 
limit. 

The motion of a pair of poles becoming attached to the imaginary axis bears some similarity to 
that found in a recent paper,'^ where as the (dynamic) spontaneous symmetry breaking happens, 
a pair of QNM poles becomes glued to the imaginary axis. In their case one of the poles stays 
at the origin, signaling a gapless Goldstone boson. We will see below one peculiar limit where a 
conductivity pole hits the origin. 



E. Truncations 



If one is interested in the behavior on the real frequency axis only, the expression for the 
conductivity arising from the AdS/CFT correspondence can be truncated to a finite number of 
poles and zeros. For instance, in a parameter regime believed to be of interest to the a wide class 
of CFTs, the conductivity has a single purely imaginary pole, accompanied by satellite poles off 
the imaginary axis. By truncating the number of poles we obtain an excellent approximation to 



the exact dependence as we show in Fig. 7(a) Up counts the number of poles/zeros, not counting 
the time-reversal partners. The truncated conductivity reads 



- n ^ (^-Cn)(^ + C) 



(37) 



where 2np — 1 is the odd number of poles or zeros (the —1 follows because the Drude pole/zero is 
its own time-reversal partner). The value of the zero is obtained by fixing a{0) = ctq. Just like 
TTo, it lies on the imaginary axis: 



Up-l 

n=l 



Cn 



(3^ 
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Exact 




FIG. 7. Conductivity (a), and its dual (b), a = 1/cj, arising from a holographic treatment with a 
truncated number of poles, 2np — 1. One pole lies on the imaginary axis, the Drude pole, while — 1 
pairs have a finite real part. The Drude form is characterized by a single pole: a = <to/(1 — itor). 



It is included so that the truncated conductivity goes to a finite constant as lim^_j.oo cr = a^o > 0. 



Fig. 7(b) shows the corresponding dual conductivity, a{w) = l/a{w), whose poles/zeros correspond 
to the zeros/poles of a. Note that the real part of the dual Drude conductivity, a = 1/a = 
(1 — w/7i^)/ao, is trivially constant (for real frequencies). 



III. EMERGENCE OF DRUDE FORM IN LARGE-iV CFT'S AND BEYOND 



In this section, we examine the conductivity of CFTs such as the critical point of the 0{N) 
model in a perturbative expansion away from the free theory obtained for = oo, with a 
focus on the emergent pole structure. We are thus approaching a general correlated CFT from 
the free quantum gas limit, as illustrated in the l.h.s. of Fig. [T| in contrast to the holographic 
approach, the r.h.s.. Our main example, though not the only one, is the NLcjM formulation of 
the 0{N) model. We show that the small-frequency quantum critical conductivity in the large- iV 
limit accurately satisfies the Drude form: 



cr W 



1 — iUT 



(39) 



The quantum Boltzmann equation (QBE) approach in the hydrodynamic regime thus captures 
the leading QNM at small frequencies, but is limited in that it misses the other poles and all the 
zeros. Although it would be desirable to have a method that captures the full analytic structure of 
the conductivity of CFTs such as the 0{N) model, the Drude pole nonetheless contains essential 
information in the d.c. limit. In addition, we can use the Drude form to verify conductivity 
sum-rules. 

The fact that a single pole can capture the small-frequency complex conductivity at large but 



22 



finite N can seem a priori surprising given that the QBE that is solved to obtain a is fairly 
complicated, including both elastic and inelastic scattering of the critical quasiparticles. Below, 
we untangle previous analyzes by providing a transparent solution to the QBE, which leads to 
the emergent Drude behavior of the low-frequency conductivity. Although we focus mainly on 
the 0{N) model, we provide similar results for a particular gauged 0{N) model as well as for a 
fermionic CFT. 



Let us first consider the case of the pure 0{N) model. We focus on the small frequency limit, 
uj where the conductivity a adopts the universal scaling fornPl 

where e is the quantum of charge, and the subscript / in the scaling function S reminds us that 
it is valid only at small frequencies, u <^T. The factors of are such that the small-frequency 
conductivity becomes a delta function at = oo, the free limit. For u <ti T, the conductivity is 
dominated by the incoherent inelastic scattering processes between the bosons. When A^ is very 
large these scattering processes can be treated perturbatively in 1/A^. Under an applied oscillatory 
field that couples to the charge, the distribution functions of the bosonic positive/negative (+/— ) 
charge excitations are modified to linear order according to f±{k, u) = nB{ek)2n6{u)+sE-kip{k, u). 
It can be shown that the linearized QBE for the deviation (p takes the formP^l 

-iojLp + g{p) = -F{p)lp + j dp'K{p,p')ip{p') , (41) 

where we have rescaled the frequency, u = Nu/T, defined the dimensionless momentum p = k/T, 
and absorbed factors of T and A^ into (f. The function g{k/T) = Td^^nsitk) acts as "source" for 
the QBE, where e\ = A(T)^ + k'^ and A oc T; we refer the reader to Appendix [A| and [b| for further 
details. Solving the equation numerically, we find that to great precision, the solution satisfies the 
simple form 

^{P,^)= (42) 
where J^{p) is a monotonous function whose behavior closely resembles that of F{p) as can be 



seen in Fig. 8(b) The case J-" = F would be the exact solution in the absence of the kernel K. 
(The latter complicates the analysis and prevents analytical solubility.) We see that the effect the 
kernel K is to renormalize F to T , which encodes all the information about the non-trivial inelastic 



scattering processes. The corresponding solution for the conductivity is shown in Fig. 8(a) it can 
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FIG. 8. a) Universal scaling function for the small-frequency conductivity, of the quantum 

critical 0{N) model. The solid lines correspond to the numerical solution of the non-trivial QBE, while 
the dashed ones to the Drude form fit. b) The momentum-dependent F{p) function entering the kernel 



of the QBE, Eq. (41), and the renormalized J- function determining the solution of the QBE, Eq. (42). 



be obtainecr^ by integrating Lp: 



oi.<^T) = '^N.'rdpt^, (43) 



where A is a momentum cutoff that is used in the numerical solution. We note that as (p decays ex- 
ponentially at large momenta, a cutoff can be safely used. Interestingly, the resulting conductivity 
is found to obey a Drude form to great accuracy: 

where S(0) = 0.085 and f = r/T = 0.775 are two universal numbers that characterize the 
entire low-frequency charge response. The former yields the d.c. conductivity while the latter is a 
dimensionless scattering rate: 

ao = ^ X iVE(O) , (45) 

^ = ^ • (46) 

The plot for the Drude form is shown with dashed lines in Fig. |8(a)[ The numerical solution and 
the Drude forms are nearly indistinguishable over the entire range < w < 14.5. The emergent 
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scattering rate 1/r gives the location of the only pole of the conductivity in this limit: 



T 

WDmdc = -i — - (47) 
I\T 

As N grows, the pole approaches the origin along the imaginary axis in the LHP; once it reaches 



it, the low-frequency conductivity becomes a delta function, as shown by the arrow in Fig. 1(b) 
The N = oo conductivity is singular and cannot be described by a meromorphic function. This 
is to be expected since it describes the transport of a free gas of bosons as opposed to a generic 
correlated CFT. 

At this point, we can compare these numerical results with those from the holographic analysis. 
In the latter, we take 7 = 1/12, which saturates the stability bound on the particle-like side 
and should be the most appropriate to compare with the almost free large- 0{N) quantum 
critical point. Indeed, the further 7 is from the bound, the closer the effective theory is to the 
strongly interacting "ideal quantum fluid" limit found at 7 = 0. At 7 = 1/12, we find that 



the Drude pole is located at WD°ude ~ — 0.26i (see Fig. 5(a) or Fig. 11(a)), which translates to 
'^Dmde = ""^^TTWp'JyjgT ^ —i3.27T. On the other hand, the Drude pole of the 0(2) model obtained 



by extending the result from the large-A^ limit, Eq. (47), is located at WDrudc ~ — 'i0.65T. The 
Drude pole from the QBE approach is thus located closer to the origin compared to the one arising 
from the holographic analysis. We thus predict that higher 1/A^ corrections to the QBE will push 
the pole further down in the LHP. This is not surprising because the extension of the large-A^ 
result to A^ = 2 yields a ratio of the d.c. to high frequency conductivities, ctq/cToo, that is larger 
than within the holographic analysis: 

^ jv=2^ 2.13 , large- AT 0{N) model (48) 

(Too L(00) 

= 1 + 47 = 1.33 , holography (49) 

N=2 

where we have used S(oo) = (1 — 8?7/3)/16 )■ 0.03998 as the large-frequency scaling function for 

the conductivity of the 0{N) model at order 1/A^, with rj oc 1/A^ being the anomalous dimension 
of the boson field.!^ It is expected that higher order 1/A^ corrections will decrease this ratio and 
will thus push the Drude pole further away from the origin. 



A. Interactions spread the weight 



Using the above quasi-exact Drude dependence, we can examine the sum rule for the low- 
frequency part of the conductivity. This is a limited version of the sum rules for the full universal 
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FIG. 9. Universal scaling function for the conductivity of the gauged 0{N) model, with damped 
gauge field. The solid lines correspond to the numerical solution of the non-trivial QBE, while the dashed 
ones to the Drude form. 



conductivity, Eq. ([T]) and Eq. The sum rule reads 

POO 

/ du^J:i{u) = nD /A = 0.1723506... , (50) 
Jo 



where we have defined the constant 

3,2 i _ ]_ 



nD= I dx[l + — \ — = 0.689403 . . . , (51) 



where = 21n[(l + a/5)/2] is twice the natural logarithm of the golden ratio. The integral 
involving the Bose-Einstein function follows simply from the expression of the conductivity in the 
free theory at = 00, see Appendix [Bj In that limit, the low frequency part of the conductivity 



reads ^ai{uj) = {TtxD /2)5{uj). On the other hand the Drude form, Eq. (44), satisfies the following 
relation: 



00 nco 



du^Ei{u)= I du^i 1 = = 0.17221... (52) 

'0 Jo [I-ITU} 2 T 

where in the last equality we have used the result given above for S(0) and f. We find that the 



emergent Drude form satisfies the sum rule Eq. (50 ) within a margin of 10 , leaving plenty of room 



for numerical uncertainty. We thus see that the interactions generated at order spread the 
weight of delta function over a finite Drude peak, whose area corresponds exactly to that of the 6- 
function of the free theory at = 00. Not only is this an excellent check on the calculation, it also 
provides a constraint between the location of the Drude pole and the value of the d.c. conductivity. 
We are effectively left with a single universal number characterizing the small-frequency behavior 
of the complex conductivity at low frequencies. 
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B. Flattening the conductivity with gauge bosons 



We now consider an interesting application of the above sum rule to a gauged 0{N) model, 
where the gauge field is Landau damped by a Fermi surface of spinons,'^^'^ which breaks confor- 
mal invariance of the critical point. This field theory was shown to be relevant to the quantum 
critical Mott transition from a metal to quantum spin liquid,'^ as well as for the quantum critical 
transition between a Neel-ordered Fermi-pocket metal and a non-FL algebraic charge liquid, called 
a "doublon metal" .'^Sl It was showrpl that the same scaling form, Eq. (40), holds as for the pure 



rotor model, Eq. (40), since only the static gauge fiuctuations contribute, the dynamical ones be- 
ing strongly quenched by the Landau damping. This phenomenon was referred to as a "fermionic 
Higgs mechanism" The numerical solution to the QBE including the static gauge fiuctuations is 
shown in Fig. M (for details, see Ref. El). As in the case of the pure 0{N) CFT, it obeys a Drude 



form, Eq. (^39^ with Drude parameters Eq. (45) and Eq. (46), this time with numerical values: 



S(0) = 0.010, f = 0.092. (53) 

The d.c. conductivity S(0) is smaller than in the un-gauged 0{N) model due to the additional 
scattering channel: the gauge bosons. The static gauge fiuctuations are actually quite strong and 
thus appreciably decrease the scattering time. The numerical solution and the Drude form agree 
very well again. Note the large range of scaled frequencies over which the agreement occurs. The 
deviations between the Drude and numerical solution seem slightly larger than in the pure rotor 
theory probably due to numerical uncertainties. The low-frequency sum rule for the conductivity. 



Eq. (l50j), yields: 

7rS(0) 



0.1720..., (54) 



differing from irD/A by only 3.5 x 10~^. We see that as we add Landau damped gauge bosons to 
the pure 0{N) model, we flatten the conductivity while keeping the emergent Drude form. The 
interactions, again, preserve the weight of the Drude peak. 



C. Fermionic CFT 

We now examine the conductivity in an interacting CFT of Dirac fermions that arises in a model 
for transitions between fractional quantum Hall and normal states.'^ The field theory consists of 
two Dirac fermions with masses Mi and M2 coupled to a Chern-Simons gauge field. The latter 
attaches fiux tubes to each Dirac fermion converting it to a Dirac anyon with statistical parameter 
(1 — a), where a = g"^ / (27r), g being the gauge coupling. The coupling a characterizes the strength 
of the long range interaction between the Dirac quasiparticles mediated by the Chern-Simons 
field. When Mi,M2 > the system is in a fractional quantum Hall state with Hall conductivity 
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FIG. 10. Universal scaling functions for the conductivity of interacting Dirac fermions a) as computed 
by solving a QBE^b) from the Drude form fitted to a). 



o'xy = e^q^/{h{l — a)), where qe is the electric charge of each Dirac quasiparticle. The transition 
to an insulating state is obtained at the point where Mi changes sign while M2 is taken to be large 
and constant. At the quantum critical point, the Mi Dirac quasiparticles coupled to the Chern- 
Simons gauge field yield a finite and universal longitudinal conductivity, whose small-frequency 



functional form is analogous to Eq. (40): 



„2 



~^fM = ^^^Z (-^ ) , (55) 



where 1/a plays the same role as N did in the 0{N) model and is taken be large. To be more 
accurate, a is the response to the total electric field, including a contribution from the emergent 
Chern-Simons field. It can be simply related to the physical conductivity.'^ The superscript "qp" 
reminds us that this is the low-frequency contribution arising from the scattering of thermally 
excited quasiparticles with each other; it is simply a different notation for E/. 

A QBE was numerically solvecPto leading order in a^, and the result is reproduced in Fig. 10 (a 



while the corresponding Drude form fit is shown in Fig. 10(b) Again, both plots agree very well. 



The two universal Drude parameters extracted from the fit are: 

E^P(O) ^ 0.437 , f ^ 0.664 . (56) 
The sum rule for the model is given in Ref. |5l 

r-3?K(^)] = ^ = 0.3466..., (57) 
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where Co = uj/^a^T). By using the Drude form E^p(ci;) = E^p(0)/(1 — iruj), we find 



r -ntlUQ)] ^ 0.33 (58) 

Jo 

The agreement is again quite good. 

In summary, we have shown that the Drude form with its single pole captures well the low- 
frequency hydrodynamic conductivity of different CFTs, a fact that was not appreciated before. 
We have also seen that such a description holds for a deformation of the 0{N) model to include 
nearly static gauge modes. Low frequency sum rules where verified in all the models, and serve as 
a useful guide in the study of interactions on the charge response. 



IV. CONCLUSIONS 

The main thesis of this paper is that charge transport of CFTs in 2+1 dimensions is most 
efficiently described by a knowledge of the poles and zeros of the conductivity in the lower-half of 
the complex frequency plane. Truncation to a small number of poles and zeros gives an accurate 
description of the crossover from the hydrodynamic physics at small frequencies to the quantum- 



critical physics at high frequencies, as was shown in Section HE Such truncated forms can be 
used as a comparison ground with experimentally or numerically measured charge response at 
conformal quantum critical points. We also showed that the conductivity of CFTs with a global 
U(l) symmetry exactly obeys two sum rules, Eq. ([T]) and Eq. (|2]), for the conductivity and its 
(S-dual) inverse. The holographic computations presented here are the first to satisfy both sum 
rules, while earlier quantum Boltzmann-theory computations satisfy only one of them. 

In the holographic approach, the poles and zeros of the conductivity are identified with quasi- 
normal modes of gauge field fiuctuations in the presence of a horizon. These quasi-normal modes are 
the proper degrees of freedom for describing quantum critical transport, replacing the role played 
by the quasiparticles in Boltzmann transport theory. We presented results for the quasi-normal 
mode frequencies in an effective holographic theory for CFTs which kept up to four derivative 
terms in a gradient expansion. 

We expect that the quasi-normal modes will help describe a wide variety of dynamical phe- 
nomena in strongly-interacting quantum systems, including those associated with deviations from 
equilibrium.^ The quasi-normal mode poles and zeros should also help in the analytic continuation 
of imaginary time data obtained from quantum Monte Carlo simulations. 
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Appendix A: Conductivity sum rules 

Conductivity sum rules are familiar in condensed matter physics in systems with a finite lattice 
cutoff. The standard derivation starting from the Kubo formula for a general Hamiltonian, "H, 
yielddSa 

poo -1 

duj^a{uj) = --\im—m,p{q)],p{-q)]), (Al) 

where p{q) is the density operator at wavevector g, and V is the system's volume. It is now our 
task to understand the structure of the commutators on the right-hand-side in the scaling limit 
appropriate for a CFT in 2+1 dimensions. 



In quantum field theory, the r.h.s. of Eq. (Al) has the structure of an ultraviolet divergent 



Schwinger contact term.^S The divergence is acceptable to us, because the sum rule in Eq. ([!]) 
is convergent only after the subtraction of the constant cxoo term. The important issue for us is 



whether the r.h.s. of Eq. ( Al ) has any finite corrections which depend upon infrared energy scales 
such as the temperature or chemical potential (/x). If such finite corrections are absent, then the 
sum rules in Eq. ([T| and Eq. ^ follow immediately, because cXoo is the value of the (j{uj) at T = 
and /i = 0, and the integral is independent of T and p. 

It is useful to analyze this issue first for a simple CFT of free Dirac fermions. Here we can 
regularize the Dirac fermions on a honeycomb lattice (as in graphene). Fortunately, such a sum 
rule analysis for the honeycomb lattice has already been carried out in Ref. |?6'. On a lattice with 
spacing a, Fermi velocity f^, temperature T, and chemical potential p, they find when T and p 
are smaller than the bandwidth that 

X = ci^ + ^/(/i/T), (A2) 
a Vp 

for some constant ci and function /. Observe that this is divergent in the continuum limit (a — )■ 
at fixed Vp, T, p), but the leading portion dependent upon T and p vanishes. So there is no 
dependence of X of the CFT upon p and T. 

Let us now carry out the corresponding analysis for the large- limit of the 0{N) rotor model. 
This is an interacting theory at finite A^, and we will see that the scaling limit has to be taken 
carefully so that we remain properly in the vicinity of the conformal fixed point in the presence of 
infrared perturbations like T or deviations from the critical point. We regularize the rotor model 
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on a square lattice of sites ij, spacing a, with the Hamiltonian 

^ = 1^ E + 17 - (A3) 

i ^ {ij) 

where (j)ia, with a = 1 . . . N are the rotor co-ordinates which obey the constraint 

J2^l = ^ (A4) 

a 

at all sites i. The their conjugate momenta with 

5 

[4>ia,T^jfi] = iSafi-^- (A5) 

The coupling constant g is used to fix the model in the vicinity of the critical point at g = gc, and 
we will take the continuum limit a — > at fixed velocity c and T. In the large N limit, the critical 
point is at 

1 C f du 1 , . , 

(A6) 



9c ifeeBZ J 27r [w^ + 2(c/a)2(2 - cos(/c^a) - cos(A;j^a))] ' 
This determines gc ~ S.llac. If we move away from the critical point, or to non-zero tempera- 



tures, then the Lagrange multiplier enforcing the constraint Eq. (A4) induces an energy gap A(T) 
determined by 

~9 ^ Lbz^? K + 2(c/a)2(2 - cos(fc,.a) - cos(V)) + A'(^)] ' ^^^^ 

where a;„ are the bosonic Matsubara frequencies. We will take the limit a — at fixed A(T) and 
T. In this limit we have 

1 1 A(0) 



9 9c 47r 



(A^ 



The density operator is 

p{q) = 5Z e'''^''''^"/? (A9) 

i 

where lap is one of the antisymmetric generators of 0{N) normalized so that Tr(/^) = — 1. Evalu- 



ating the commutator in Eq. (Al) we find 



m, p{q)] , p{-q)] = -—Y. - e*'"" I' • (MO) 



9 
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So taking the limit the long wavelength limit yields 



lini — [[n, p{q)] , p{-q)] 



(All) 



Using Eq. (A4), we can now write the conductivity sum rule as 



I 



2g AgV 

TTC^ Hi? 



fceBZ 



(2 — co^ikxo) — cos^kyo)) 



[ul + 2(c/a)2(2 - cos(fc^a) - cos(A;ya)) + A^] ' 



(A12) 



Evaluating the frequency summation, and then taking the limit a — > 0, we obtain the expansion 



TIC 

2^ 



tti — V CX2 — a — a — - 



oo t2 



+ . . . 



(A13) 



where ai ~ 0.75 and a2 ~ 0.13. The crucial feature of this result is that there is no term ~ A, 
and all terms containing A vanish as a — )■ 0. A term ~ A does appear if we choose a general A 



which does not obey Eq. (A7) and then evaluate Eq. (All). Thus the imposition of the constraint 



Eq. (A4) at all T was important for the absence of such a term. The general features of Eq. (A13) 



are similar to Eq. (A2), and so the same conclusions apply. 



Appendix B: Analytic structure in the N ^ oo limit of the 0{N) model 

This appendix notes a few features of the conductivity of the 0{N) rotor model in the complex 



frequency plane, in the N oo limit. For the model in Eq. (A3), the conductivity as a function 
of the complex frequency z follows from Ref. H) 



iTD 



iz 
An 



Q?(z^ 



(Bl) 



where the contour of Vl integration determines the specific choice of the current correlator, and the 
Drude weight scales linearly with the temperature. We have defined the numerical constant 



D 



fisinh2(fi/(2T))' 



(B2) 



whose value is given in Eq. (51) 



The retarded response function a^{z) is obtained by choosing z in the UHP, and the contour 
of integration along the real frequency axis. This function a^{z) is analytic in the UHP, and 
has a pole at z = and branch points at z = ±2A. We can perform the analytic continuation 
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of a^{z) into the lower-half plane by deforming the contour of Q integration into the lower-half 
plane, so that it is always below the points ±z/2. Because of the presence of these branch points, 
the analytic continuation of (t^{z) into the lower-half plane is not unique, and depends upon the 
path of z around the branch points. This is a key difference from the holographic results of the 
present paper, which had no branch points and a unique analytic continuation into the LHP. We 
expect that fully incorporating corrections will make the 0{N) model result similar to the 
holographic computation. We have already demonstrated this for the case of the pole at 2; = 0, 
which becomes a LHP Drude pole. However a careful analysis of corrections determining the 
fate of the branch points a.t z = ±2 A has not yet been carried out. 



In any case, the physical value on the real axis cx^lu + iO^) is unique, and was shown in 
Fig. 3(d) At the critical point, this is to be evaluated at A = 0T, where 6 = 21n((v^ + l)/2). 
Curiously, for this value of A, we find zeros of the conductivity on the real axis branch points, 
with a^{±2QT + iO^) = 0. So the structure of poles and zeros of the N = 00 conductivity has 
a remarkable similarity to the 7 > holographic results, as was reviewed in Fig. |3j The pole at 
z = of the N = 00 theory corresponds to the closest pole on the negative imaginary axis of the 
holographic result, as we have already noted. And the zeros ai z = ±20T of theory correspond 



to the two zeros closest to the real axis in Fig. 4(b) 



Finally, we can verify that the sum rule in Eq. ([T]) is satisfied by Eq. (Bl) 

1 



3f?(T^(u; + iO+) - 



16 



0. 



(B3) 



where we have used (Too = 1/16. Note that this result is obeyed only for A = OT, and not for 
other values of A, as is expected from the considerations in Appendix |Aj Also, as noted in the 
introduction, the inverse sum rule in Eq. ^ is not satisfied by Eq. (Bl); indeed, in this case the 
integral diverges because of the real-axis zeros of the conductivity. 



Let us also mention that the analytic structure of response functions of the 0{N) model was 
also examined recently in Ref. ITflaway from the CFT critical point, but at T = 0. In the ordered 
phase with broken 0{N) symmetry, poles were found in the lower-half plane corresponding to the 
Higgs excitations damped by multiple spin-wave emission. 
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Appendix C: Differential equation for the numerical solution of the conductivity 



We first factor out tlie singular part of Ay near the horizon: Ay = {1 — u) ^^F{u). Making this 



substitution in the EoM for Ay, Eq. (12), we obtain the following differential equation for F: 



= F"- 



3^2(1 -4(1 - 2u^)-f) 



2iw 



F' 



+ 



(1 -n3)(l + 4^37) l-u^ 
iw {{1+U + u^){l + 2u + 4^2(3 + Au + 5u^)-f) -i{2 + u){A + u + u'^){l + Au^'y)w) 

(1 - m)(1 + M + m2)2(i + 4^3^) 



(CI) 



This is to be compared with the simpler form of the equation for the full Ay, Eq. (12). The two 
boundary conditions at the horizon read 



m 

F'{1) 



1 



iw{i + 2w + 87(2^ + w)) 
(1 +47)(i + 2w) 



(C2) 
(C3) 



The second condition follows from the solution of the differential equation near u = 1: F{u) 



1 — (1 — u)r, with r being the r.h.s. of Eq. (C3). The numerical solution is shown in Fig. H^and 



in Fig. ITTl where the poles and zeros in the LHP can be seen more precisely. 



Appendix D: WKB analysis for asymptotic quasi- normal modes 

The goal of the WKB analysis is to identify the QNMs of the gauge field at large frequencies, 
\w\ ^ 1. According to the AdS/CFT correspondence, these frequencies can then be put in corre- 
spondence with the poles of the gauge correlation function Qyy proportional to the conductivity. 



Eq. (14). The standard analysis examines the solutions to Eq. (D15) near 1) the black hole sin- 
gularity, 2) the event horizon, and 3) the asymptotic boundary. Matching of the solutions usually 
gives an expression for a set of discrete QNM frequencies. Generically one obtains two solution for 
Ay, with one vanishing as the boundary is approached. Discarding the non- vanishing one leads to 
a "quantization" condition on the QNMs. 

As mentioned in the main text, the EoM for the y-component of the gauge field reads: 



h' q?;;2 



h f 9 

The second equality follows from h = fg. We can change coordinates to bring this equation into a 
Schrodinger form, which will be more convenient for the analysis of the QNMs. To do so, we want 
to transform away the linear-derivative term. One way involves changing variables to dx = du/f. 
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(a)5i{a(7 - 1/12)} (b)5i{a(7 = 1/12)} 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

w' w' 



(c)5R{a(7 = -1/12)} (d)K{a(7 = -1/12)} 

FIG. 11. Conductivity a and its dual a = \/a m the LHP, w" = < 0, for I7I = 1/12. There is a 
quahtative correspondence of the pole structure between a{'w;j) and a{w, —7). Note that the poles of 
a{'W]'y) are the zeros of a{w,'y). 



as we illustrate below. 

Before going into the WKB analysis, let us first review the simplest scenario, 7 = 0, i.e. in the 
absence of the function g arising from the Weyl curvature coupling. The exact solution is obtained 
by using the new (complex) coordinate z: 

dz 3 3 , 



This puts Eq. (Dl) in the form: 

d^Ay + w^Ay = (D4) 
with solutions: e^™^. To apply the boundary condition we need to examine the explicit form of 
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z(u). Integrating Eq. (D3), we obtain 



3 

z{u) = Y.j;—^Hl-u/up) (D5) 

p=l ^ ^'^p' 

where Up are the 3 zeros of /. They are simply the cubic roots of unity: = 1, i.e. 

Ml = 1 (D6) 
U2 = -{l + tV3)/2, U3 = u; (D7) 

which is trivially found by noting that 1 — -u^ = (1 — u){u^ + « + !). We give a few properties 
of the generating polynomial / and its roots that will be useful for future analysis. First, the 
derivative of / permutes U2 and while leaving ui invariant (up to signs): /'(mi)/3 = —Ui and 
/'(u2)/3 = — M3. As a result, we get the following identities: 

3 

^ Mp = (D8) 

p=i f'^'^p) I otherwise 
Recall that we need to apply an infalling boundary condition. Ay ^ (1 — near the event 



horizon, u = 1. Using Eq. (D5), we find that as m — )■ 1, 

^±iwz ^ X (1 - u)^'"" (DIO) 

where C± = e^^'^{^'^^+'^/^)/'^_ Hence, the boundary condition selects Ay = e™^ . This in turn yields: 
a = ~^^^|„_5.o — ~^ 3wii-u3) L^o ~ ^- expected the conductivity of the CFT holographically 
dual to the Einstein-Maxwell theory on S-AdS4 is constant for all complex frequencies, hence self- 
dual. We now include a finite 7, which prevents analytical solubility, just like the collision 
term did for the 0{N) model. 



We wish to transform Eq. (Dl) into a Schrodinger form. To facilitate comparison with the 
literature, notably with Ref. US] which serves as a guide for our analysis, we shall perform the 
WKB analysis starting with the coordinate r = 1/u instead of u. This is the radial holographic 
coordinate introduced in the main body, with the difference that it is rescaled by tq. We define 
f = r^/ = — r^^, and the corresponding new tortoise coordinate (the analogue of z introduced 
above): 
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In terms of x, the EoM for Ay becomes: 

(fAy 1 dq dA o , ^ „ 

--r^ + --r^ + ^^v = ^^ z/ = 3«7. pl2 
dx"^ g dx dx 

We have defined the rescaled frequency v to simphfy the comparison with previous works. We 
note that in the limit where 7 = 0, the hnear derivative term vanishes and we are left with a trivial 
harmonic equation as above. For finite 7, we can remove such a term by introducing two functions 
to parameterize Ay-. 

Ay = G{x)%l){x) , (D13) 

where in order for if) to satisfy an equation of the Schrodinger form, G needs to satisfy the first 
order differential equation: 

dG 1 dq ^ s 
dx 2g dx 



This can be solved in general by G = l/^/g = l/yl + A^yu^. The resulting "Schrodinger" equation 
for is: 



d'^ 



dx"^ 

where 



+ W{x)i) = u^i), (D15) 



^ = ;^p^[2-' + - 5)r^ - 147] . (D16) 
The potential W prevents the exact solubility of the equation, and as expected vanishes as 7 — > 0. 



In that limit, G — )■ 1 and — )■ 0, and the equation reduces to the harmonic one Eq. (D4). Note 
that the potential vanishes at the boundary, r = 00, just as the Weyl curvature does. 



The underlying idea of the WKB method is to examine the behavior of Ay or on the Stokes 
line in the complex r-plane defined via: 

S5(z/x) = 0. (D17) 
The first step is thus to identify this Stokes line by studying the behavior of the tortoise in terms 



37 




FIG. 12. The Stokes line, = 0, in black in the complex r-plane; r = corresponds to the 

intersection point of the two branches of the Stokes line. The color shading represents the value of '^{ux). 
The 3 branch cuts coming from the logarithms are clearly visible. 



of r. As above, the defining relation for the tortoise can be integrated to give: 

1 ^ 1 



p=i 



= - [ln(l — r) + a* ln(l — a*r) + a ln(l — ar)] 
3 

where ri = l,r2 = a,r3 = a* = are the three cubic roots of unity, with a 
precisely the Up introduced above. Near r = 0, oo, the tortoise scales like 



X ~ — - 



r -> 



X ^ Xq , r— i-oo 

r 



(D18) 
(D19) 

(D20) 

(D21) 
(D22) 



respectively, where we have introduced 



Xq = x{r — )■ oo) 



( 

9 



(D23) 



which will play a central role in the WKB analysis.!^ Its value is well-defined due to the absence 
of monodromy at infinity, even in the presence of the 3 branch cuts coming from the logarithms, 
see Fig. 12 The value of xq dictates that of u via vxq G M: = Ce*'^/^, where C ^ I^- Iii particular, 
from this and Eq. (D21), we see that the branch of the Stokes line that extends to infinity follows 
the line r = pe"^/'^, where p is real. Near the origin, we have '^{e^'^^^x) ; 



-3(e'^/V2)/2, which 

p G M, in addition to r = pe*'^/^. These two branches of the Stokes line cross at 



implies r = pe 

the origin as we show in Fig. 12 We now proceed to the WKB analysis by examining the solution 
to Eq. (D15) in the vicinity of r = oo, 0, 1. 
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Near r = oo: The potential W{r) is irrelevant since W ^ 1/r. This is not surprising since we 
expect 7 to be irrelevant near the UV boundary and W oc j. The equation becomes harmonic. 
We write the solution in terms of the shifted variable, x — xq, and use Bessel functions although 
simple sines and cosines would suffice; this allows us to compare with other QNM analyses.'^ We 
have 



ip{x) = B+^/2nu{x - XQ)Jj^/2{y{x - xq)) + 5_ v^27n^^ar^xJ^J_j^/2(z/(x - xq)) (D24) 



where j^o = 1, and Ji/2{z) = ^/2J^: mi{z) / y/z , J-i/2{z) = A/27vrcos(2;)/A/^. As we have discussed 
in the main text, we need to impose the vanishing of Ay = tpG at the boundary, which leads to 
ip{xo) = since G{xo) = 1. We thus have our first constraint, = 0. 
Near r = 0: Near the black hole singularity, the potential diverges: 

with jo = 5/2. In the second inequality we have used x ~ — r^/2 near the singularity. We thus 
have the Bessel solution: 

iIj{x) = V 27rz/x Jjo/2 (z^x) + A_\/2jnJxJ^jg/2{i'x) (D26) 

We can match the solutions near r = oo and r = using the asymptotic expansion for z ^ 1: 
Ja{z) ~ ^y2/{7rz) cos[z — (1 + 2a)7r/4]. Extending starting near the origin r = 0, 

iIj{x) ^ 2A+ cos(z/x - a+) + 2A^ cos(z/x - a_) (D27) 
= (A+e-*"+ + A^e-'""- + (A+e*"+ + A_e*"-)e-^^^ (D28) 

where we have defined a± = {1 ± jo)^/4. On the other hand extending from r = oo we get 

ip ^ 2B+ cos[u{x - xo) - (D29) 
= 5_^e-*^+e*''(^-^») + 5^e*'^+e-*''(^-^°) (D30) 

where /3+ = 7r/2. Matching both solutions by equating the ratios of the coefficients of e^*^^ yields 
another constraint: 

A^ sin(z/xo + /3+ — «+) + A^ sin(z/xo + (3^ — «-) = . (D31) 

Near r = 1: We then want to match the behavior on the Stokes branch r = pe™!"^ with that 
near the black hole event horizon r = 1. First, we have the small-2; expansion Ja{z) ~ z°'w(^z\ 
where w{z) is an even and holomorphic function, w{z) = o-Fi(a + l; — z^/4)/(2'^r(a + l)), where 0-^1 
is an instance of the hypergeometric function. We will rotate from the branch r = pe*'^/^, p E 
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to r = pe ^ p G M"*". Using a; ~ near r = 0, the 7r/2 r-rotation becomes a vr x-rotation: 



V27re-^^i/xJ±j„/2(e"*^z/x) = e-^(i±Jo)^/V27rz/xJ±jo/2(«^a;) (D32) 

^ 2e-*2°± cos(z/x - a±) (D33) 

Using this we have the following behavior on the r = pe~*'^/^, p G branch: 

~ 2A+e-*2°+ cos(-z/x - a+) + 2A„e~*2"- cos{-ux - a„) (D34) 

= (A+e-*"+ + A_e-*"- )e^^^ + (^+6"*^"+ + A.e"*^"- )e-^'^^ (D35) 

We know that at the horizon, ip^x) ~ e*'^^ in order to satisfy the infalling condition, consequently 

A+e-^^°+ + A„e-^^°- = . (D36) 



Combining Eq. (D31) and Eq. (D36), we find get a condition that the homogeneous system of 



equations needs to satisfy in order to have a solution: 



-i3o_ 



det 



sin(z/xo + /3+ — a+) sin(z/xo + /?-( 



(D37) 



This equation leads to the general solution for the asymptotic QNMs: 

3w7Xo = ^ - 27rn , n G N & n > 1 
where we have switched back to w = z//3. We find two solutions for the offset parameter ^: 



(D38) 



^1 = 2z tanh 



v^+(-l 



^2 = 2 tan"^ ^ ^ '- ' ^ 



-2.356 -zO.173, 



0.785 - zO.173 



(D39) 
(D40) 



The offset and gap, defined via w = [gap] — n [offset] for large n, are given by 



offset 



gap 



3x0 ' 
27r _ 

3xn 



(D41) 
(D42) 



where the offset obtained using ^i_2 is — 0.283 — iO. 586 or 0.150 + iO. 164, respectively. Interestingly, 
we note that these results for the asymptotic QNMs are independent of the value of 7, as long as 



it is finite. In contrast, if 7 = 0, we obtain jo = joo = 1, and the determinant condition Eq. (D37) 
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leads to e*^^° = 0, which has no finite solution. This is in agreement with the exact solution: there 
are no QNMs when 7 = because the corresponding conductivity is a constant function. 
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